clc;
phi=pi/3;
T=1e-3;
Td=1/44.2e6;
t=Td:Td:T;
Omega0=2*pi*1.5e6;
Omega=2*pi*[-3000:10:3000];
figure (1)
plot (cos(Phi))
A=5;
NPSP=10; % число битов модуляции на мс
TPSP=1e-3;% длительность соответ NPSP
h=sign(randn(1,NPSP)); %зададим линейнонарастающую фазу
figure (2) 
plot (h)
dt=0;
ind=fix(mod((t+dt)/TPSP,1).*NPSP)+1;%округляем в болшую сторону,строим ступеньку
hop=sign(randn(1,NPSP));
dtop=0;
indop=fix(mod((t+dtop)/TPSP,1).*NPSP)+1;
Omegaop=2*pi*5;
Phiop=Omega0*t+Omegaop*t;
figure (3) 
plot (Sop)

Sop=hop(indop).*cos(Phiop);
SopQ=h(indop).*sin(Phiop);
for i=1:length(Omega)
Phi=Omega0*t+Omega(i)*t+phi;
S=A*h(ind).*cos(Phi);
I(i)=S*Sop.';%получили синфазную и квадратурную компоненту
Q(i)=S*SopQ.';
K(i)=sqrt((I(i))^2+(Q(i))^2);
end

figure (4) 
plot (I)
figure (5) 
plot (Q)
figure (6)
plot (Omega,K)